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I. INTRODUCTION 

The investigation of the dynamics of coupling between 
lower and upper atmosphere through wave generation, wave 
propagation and wave dissipation is essential for the under- 
standing of the motion of the atmosphere, in particular, 
accurate information about atmospheric disturbances and wave 
propagation is required for the development of a model of 
the perturbed neutral atmosphere for the use in orbital/ 
environment mission performance assessment of the Space Shuttle 
vehicle . 

Coupling between ffe lower atmosphere and the ionosphere 
can be detected from ionospheric wave -like disturbances which 
are frequently observed on ground-based ionospheric sounding 
records as perturbations in the electron densities. Recently, 
Georges (1968, 1972), Baker and Davies (1969), and Davies and 
Jones (1972) reported observations of ionospheric wave-like 
disturbances which showed some correlation with severe weather 
activity, and suggested that infrasonic waves could be generated 
by the storm systems. Hven though severe weather activity as 
possible sources of gravity-acoustic waves has been suggested, 
additional and more confirmative experimental evidence to support 
such a hypothesis is urgently rocessary. 

To provide more experimental evidence on the dynamics 
of upper atmospheric disturbances during '•evere weather activity, 
a continuous wave-spectrum, high fretiuency, Dopi'ler sounder 
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array has been employed. The purpose uf the study is to 
investigate the coupling between tne lower and upper atmo- 
sphere (ionosphere) during severe weather activity. 

Experimental observations of ionospheric wave-like 
disturbances during severe weather activity are from the NASA/ 
Marshall Space Flight Center CW Doppler sounder array located 
in the North Alabama region. Since January 1973, we have 
continuou;: ly monitored tlie Doppler frequency fluctuations of 
nine field transmitters nt three stations which cover an area 
of 150 Km^ 

In Chapter II, wc describe the cliaracterist ics of 
atmospheric waves. In (Jiaptcr HI, experimental operations 
of Doppler sounder array which h-ive been used to detect 
ionospheric disturbances are discussed. The procedures of 
data analysis from the rcTO’^ds of Doppler fluct’.iations are 
presented in Chapter TV. In Chapter V, five events of severe 
thunderstorms and tornadoes, including the extreme tornado 
outbreak of April 3, 1974, were iinaiyzed. In Chapter VI, we 
have discussed the effect of Doppler shift due to eastward 
wind which may affect the v;ave-period of acoustic-gravity waves 
observed during severe weather activity. Furthermore, the 
criterion for the excitation of S- shaped record or Doppler 
"fold-backs," which is very often observed <Diring tornado 
activity, is derived and is compared with possible tornado 
conditions . 
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II. ATMOSPHERIC WAVES 

Typical wave motions in the neutral atmosphere are the 
propagation of sonic waves. The normal frequency range of 
the sonic spectrum of the human ear is 20 Hz to 20 KHz. Based 
on this frequency range, the sonic waves are generally divided 
into the following three groups: 

(1) Audio-Sonic: Waves 

The frequency range is within the spectrum of 
human car. 

( 2 ) Ultrasonic Wave s 

The frequency range is above audio-sonic waves. 

The upper frequency limits for ultrasonic waves 
approach the relaxation frequencies of metals, around 
10® MHz, beyond which materials can no longer respond 
to the input of mechanical wave energy. 

(3) Infrasonic Waves 

The frequency range is below audio-sonic waves. 

In the real atmosphere, the presence of gravitational force 
causes the density to be vertically nonuniformly distributed 
which is responsible for the buoyancy oscillations of the fluid 
parcel of the atmosphere. The buoyancy frequency, which is also 
termed Brunt -Vaisala frequency, is given by 
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“B ' 'f fd * (2-1) 

where g is the gravitational acceleration; T, the temperature, 
z, the spatial coordinate along the vertical height; and r^, 
the dry adiabatic lapse rate which is defined 

Td ' I (^) . (2-2) 

Here R is the universal gas constant and y * 
is the specific heat at constant pressure, and is the 

specific heat. The height variations of molecular temperature 
T, sound speed a, ratio of specific heat y, and Brunt-Vaisala 
wave period tg (= Zw/wg) arc given in Figures 1, 2, 3 and 4, 
respectively. It can be seen from Figure 4 that the Brunt- 
Vaisalii frequency is practically less than 0.05 Hz in the 
atmosphere and it is therefore more convenient to speak of the 
period of the wave. This characteristic frequency caused by 
the buoyancy force is so small, as the period is several minutes 
to more than ten minutes in real atmosphere. Thus, only 
comparable low frequency waves would be affected by the 
buoyancy force. It is for this reason that the waves modified 
by the buoyancy force are limited to infrasonic waves in the 
atmosphere. In fact, the modification is so drastic that the 
resulting waves arc called acoiistic-gravity waves. 

The characteristics of acoustic-gravity waves can be 
shown in simple calculations. To do so, the atmosphere is 
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“Tempsrahure , 

Figure 1. The height vari ition of molecular temperature in 
real atmosphere. 
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Sound Speed ^ m/sec 


Figure Z. The height variation of sound speed in real 
atmosphere . 
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0 3 4- 5 6 7 8 9 /O II /2 13 14 

Brant -‘Vai5&l& ( or Buoyancy) 

Wave Periods , min. 

Figure 4. The hci;^ln vari.ition of huoysney or Brunt -VLisSlS 
wa\ •• -per «u.l in real atmosphere. 
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assumed to be nonrotating (for a wave period shorter than 
two hours) with no large-scale neutral wind. The equations 
are the equation of continuity , the equation of motion, and 
the equation of energy which are replaced by the adiabatic 
equation in this case, namely, 

|r + H • (pv) ' 0 C^*3) 


P 



2P ♦ Pg 


(2-4) 


DP 

Ut 


5t 


(2-5) 


where the convective derivative is D/Dt • 3/3t V • 7. Here 
p denotes mass density; V, velocity; P, pressure; a » (3P/3p)*/*, 
sonic velocity; and g * (0, 0, -g) , gravitational acceleration. 
For a static, motionless, isothermal atmosphere the mass density 
must be distributed exponentially as a function of height 

<p> - p^ e'*^“ (2-6) 

where p^ is the density at z » 0, and H = u^/yg (y is the ratio 
of constant pressure specific heat to constant volume specific 
heat). The corresponding pressure distribution is 

<p> « e'*/“ . (2-7) 

Let the isothermal atmosphere be perturbed according to the 
following scheme: 
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C2-8) 


p(r, t) 


> 

<p(7.)> ♦ <5o(r. t) 

P(r, t) 

— 

<P(z)> ♦ 6P(r, t) 

V(r, t) 

^ / 


6V(r, t) 


where < > denotes the equilibrium quantities and 6 implies the 
perturbations. After lin'jarizing the equations, and assuming 
that the perturbation is proportional to 


exp i (dJt - 
we have the following 


k • r) 

dispersion relation 


k* 

X 


-► 


Cl-co^V«*>*)/(l-o)j^Va)2) 


1 • u) Va;* 

tl 


(2-9) 


where k ■ w/a; ui = YS/2a = a/2H =* acoustic cutoff frequency; 

O cl 

1/2 

and ti)^ ■ (yl) P,/'i ■ buoyancy frequency at isothermal atmosphere, 
namely, 


0 ), 


U) 


BldT 

ai 


( 2 - 10 ) 


Let us discuss the properties of the dispersion relation 
(2‘9) in the following three different ranges: 


(1) li) < U)^ 

This is gravity wave. The phase velocity of gravity 
waves is always less than the velocity •: lund 

waves . 

(2) u). <(jj<w 

D P 

This is cutoff region where no acoustic-gravity 
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waves exist under the condition of isothermal 
atmosphere . 

( 3 ) 03 > 

a 

This is acoustic wave. The j lasc velocity of 
acoustic waves is equal or greater than the 
velocity of founil waves. 

In the real atmosphere, the assumption of isothermal 
atmosphere is no longer valid. This implies that the cutoff 
region of the form discussed for isothermal atmosphere no 
longer exists. Therefore, tljc overlap between the regions of 
acoustic waves .and gravity wnv< •- v/il. I be expected. 
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EXPHRIMliNTAL OPERATIONS 


The continuous wave -spectra Doppler sounder array was 
developed by Davies fl962). Davies and Baker (1966) gave a 
description of the theory and an interpretation of the data. 
NASA/Marshall Space Might Center CM Doppler sounder array 
consists of three sites with nine field transmitters operating 
at 4. 0125, 4.759 and 5.754 MHz with receivers located at 
NASA/Marshall Space Flight Center (NASA/MSFC) , Alabama (see 
Figure 5 for map description). The geographic locations of 
the field transmitters are shov/n in Tabic 1. The received 
signals from all nine transmitters are mixed with stable 
oscillators of the same nomin.il frequencies, but are offset 
several Hz higher. The heat frequencies are then recorded 
on a slowly moving (one inch per minute) magnetic tape. For 
data analysis, the tapes are played back at higher speeds 
into an audio spectrum .iiialyzer. The output of the spectrum 
analyzer is connected to a facsimile recorder which gives a 
strip chart record of the Doppler shifted frequency components 
against time. 

A block diagram of the equipment for one transmitting 
station is shown in Figure 6. Th^^ee ovenized crystal 
oscillators arc built into one chassis, as shown in the 
block diagram of Figure 7. Since each of th'' three transmitters 
must be slightly different so they can be distinguished from 
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Figure S. Locations ol’ thioc sites of field transmitters 
and receivers of ^'ASA/^far?hall Space Flight 
Center Continuous wave spectrum, high 
frequency Doppler sounder array system. 
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TABLE 1 

Geographic Location of the CW Doppler Sounder Array 


Station 

Latitude 

Longitude 

NASA/MSFC, Alabama 

34® 39'N 

86® 40'W 

Ft. McClellan, Alabama 

33® 44'N 

85® 48'W 

TVA Muscle Shoals, Alabama 

34® 46'N 

85® 38’W 

TVA Nickajack Dam, Tennessee 

35® Ol'N 

87® 38'W 
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Figure 6. JUock 'If if.ran of rranrimlnLirig 
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Block diagram of Oscillator 
chaala. 
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each other at the receiver. The actual frequencies for the 
field transmitters have been set as in Table 2. At the receivinj', 
station, the oscillators have been set exactly at the nominal 
frequencies so that in the absence of Doppler variations, the 
offsets of the transmitter oscillators are equal to the beat 
frequencies . 

The transmitters are modified versions of the Heathkit 
DK 60B. The transmitting antennas are halfwave dipoles. In 
order to maximize vertical radiation and minimize horizontal 
radiation, the height of the dipoles above ground should be 
between one-eighth to one-fourth wavelength. .Since all three 
antenna at a site are in one array, a height of approximately 
30 feet has proven to l)c practical. To null the horizontal 
radiation in the direction of the receiver, the transmitting 
antenna at the sites are oriented so that the array is broad- 
side to the receiving site and the ends of the antenna point 
to the receiver site (NASA/MSPC) . Ground waves are not expected 
to be a problem since the sites arc farther than 40 Km from 
the receivers. If a ground wave is received, it can be 
identified by its stable beat frequency on the Doppler records. 

A block diagram of the receiving and recording system is 
shown in Figure 8. The reflected ionospheric echoes from 
the nine field transmitters arc received by three Hammarlund 
SP-()00-.JX- 1 7 communication receivers at three different nominal 
frcquenc i'’s . TIics'* receivers use single conversion for channel 
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5.734 MHz 
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r.ii'i +if)iiz 
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frequencies below 7.4 MHz, and double conversion above 7.4 MHz. 
In either case, the internal crystal HF oscillator is used 
instead of the variable frequency oscillator. It is not 
necessary that the HF oscillator in the signs'*, receiver has 
the stability of the reference receiver to be mixed with the 
signal from the reference oscillator. Thus, any variations in 
frequency disappear when the signal intermediate frequency (IF) 
is mixed with the reference IF. 

In order to obtain the accurate timing of ionospheric 
disturbances, a time code is recorded on channel 4 of the data 
tape. When the tape is played back for spectrum analysis, the 
channel 4 output is mixed with whichever data channel is being 
played back, so that data and timing appear simultaneously on 
the records. In other words, channels 1, 2 and 3 are used for 
data, with a timing signal recorded on channel 4. The data 
tapes are played back on a four channel tape recorder at a 
speed of 15 inches per second. Each data tape is played back 
three times, with each playback, one data channel is mixed 
with the timing channel. The channels are recorded as follows: 
Channel 1: 4.0125 MHz, low frequency 

Channel 2: 4.759 MHz, middle frequency 

Channel 3: 5.734 MHz high frequency 

Channel 4: Timing signal. 

The output from the tape recorder is spectrum analyzed 
using a Federal Scientific UA-6B, Ubiguitous audio spectrum 
analyzer. The output from the spectrum analyzer is fed to a 
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Gift GPR facsimile recorder which j’ives a record of frequency 
against time. 

In the present study, we are intcic'^ted in the analysis 
of atmospheric waves generated by severe weather activity. A 
few examples during severe weather activity, pnrticularly 
thunderstorms and tornadoes, have been chosen ftr data analysis. 
The procedures of data analysis will bo discussed in the following 
sections. 
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IV. PROCEDURES OF DATA ANALYSIS 

The visible record of the Doppler data is obtained by 
frequency analysis of the signal. The frequency analysis is 
carried out by using a UA6B Ubiguitous audio spectrum analyzer. 

The output from the spectrum analyzer is displaced on a facsimile 
chart. In the present study, we are particularly interested in 
the investigation of the specific atmospheric- ionospheric cou- 
pling during thunderstorm and tornado activity which can be 
identified on the Doppler record. After visual identification 
of the events due to severe storm activity, the segments of the 
data tape containing the events are played back at a slower 
speed for better time resolution. By playing back the tapes 
at a slower speed and recording the data, one can achieve a 
time resolution of about five seconds, depending on the time- 
scale of interest, to study the fine structure in the Doppler 
variation. 

The Doppler fluctuations printed in facsimile chart by 
playing back the tapes arc then subjected to numerical digitizing. 
The digitized data will be used to determine time displacement 
and wave frequency of the fluctuations. In other words, for 
a given sample of the Doppler sounder record, the time displace- 
ments are determined by comptiting the cross -cor rcl at ion 
functions. Consider two time s'*ries functions given by ft) 
and XoftJ. The cross-correlation function is given by 

P 
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■ tV ” ^ I -'Tldt f4.1) 

-T 

where t is the time lag, or the time ; , ' olacemcnt between the 
two functiciis being correlated. In r; . pecial case when 
Equation (4.1) boi omes the -'Uix' ^ ^ lation function 

-T' 

<^<.0 *t'." a ( C4.2) 

J ~x 

• I 

In general, the correlation matrix [C] can be expressed by 

Clift) 

[C] = C2^(x) C^^Ct) CjjCO 

C^lft) 

In our case, if M, N and F refer to the ionospheric 
rvlfection points to the transmitters located at Muscle 
Shoals, Alabama; Nickajack Dam, Tennessee; and Port McClellan, 
Alabama to the receivers at NASA/MSFC, respectively, the 
correlation Matrix becomos 

(Cl - C^„fT) C^J,(T) 

Cp^C^) Cpp(i) 

The power spectral density anil cross power spectral 
density arc obtained by Fourier transforming the auto and 
cross correlation functions, respectively, i.e.. 
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( 4 . 5 ) 


( 4 . 6 ) 


where v ® 2nf is the angular frequency in radian/sec and f is 
frequency in Hz. In general, the cross-power spectral density 
is complex, 


where P^g(u)) is the real part and P^g is the imaginary part of 
the cross spectrum function. The phase density spectrum can be 
obtained from 


0(<*)) = 6g(u) - 0^(‘*>) 


» tan 


- 1 


’’aeW 




«. 8 ) 


where 6(j(‘*)) is phase of time series x^(t), and 6g(a)) is phase 
of time series Xg(u>), The magnitude of the complex continuous 
spectrum, called the amplitude density spectrum, is 


e 


■1 


‘/2 


( 4 . 9 ) 


The power can be calculated from the integration of power 
spectral density 
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The coherence estimates are given by 


Y(w) 




(P (u) 

^ ant 




(4.11) 


where y(w) is the complex coherence. The amplitude of the cross 
spectrum at each frequency is the product of the corresponding 
amplitudes in the two time scries. If a particular spectral 
component is absent from either time series, then it will be 
absent in the cross-spectrum and the coherence will be minimum 
at that frequency. 

In our case the frequency of the ionospheric disturbances 
caused by severe wcatb(r activity can be computed from both 
power spectral density and cross spectral density. The time 
delay between the reflection points can be obtained from cross- 
correlation function. Let and be the time delays 
between the reflection points M-N and M-F, respectively; and 
and Dj^p be the distance between the similar points, respec- 
tively. The appnrant trace velocities and Vj^^p are given 
by the following relations (See Figure 9): 


D, 


'MN 


MN 

^MN 


D, 


and 


'MF 


MF 

^MF 


(4.12) 


The magnitude horizontal trace velocity and the direction 
(|) (measured cast from north) can be determined from the following 
set of equations: 
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Figure 9. Graphical representation of wave front moving 
across the doppler sounder array. 
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rnt t, = ^MN '^'MN ' ^MF ’^MF 

"^IN '*'mN ■ '^MF MF 

and the magnitude of V^ is given by 

V* = ^MN ^MF ^'^MN ' *^MF^ 

'^MN‘"'^MF ■ ^ ^MN '^MF ^’^MN'^MF^ 


(4.13) 


(4.14) 


where and and the angles between the geographic north 

and lines MN and MF, respectively. The true velocity may be 
determined by computing the vertical trace velocity V^. The 
vertical trace velocity may be determined from the time 
displacements at various frequencies and the difference of 
reflection points for a given transmitter and receiver location, 
are assumed to lie on a vertical plane. The heights of reflections 
for various frequencies involved are Hctcrmined from the ionogram 
analysis. From the values of horizontal phase velocity V^j and 
vertical trace velocity V^, the true phase velocity may be 
determined as shown in Figure 10. 

Fortran listing of comiiutcr program for the computation 
of cross correlation functions, auto correlation functions, 
power spectral density, cross power spectral density, cross 
spectral amplitudes, cross spectral phases, cross spectral 
coherences, etc., arc included in the appendix. 
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Flgur«10. Graphical rapreaentatlon for computing phase velocity 
of Che wave front. 
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V. GHNI-RATION OF ACOUSTIC GRAVITY WAVES BY 
SEVERE WEATHER ACTIVITY 


In this section, we present the experimental evidence 
of wavelike disturbances in the ionosphere caused hy severe 
weather activity, particularly due to thunderstorms and 
tornadoes. The atmospheric waves observed in the ionospheric 
height are in the category of acoustic-gravity waves. As 
described in Section III and Section IV, the data obtained 
from the Doppler sounder are then subjected to power spectrum 
and cross power spectrum analyses. The results of experimencal 
evidence is shown as follows: 

A . Thunderstorm Observati on 

Two thunderstorm events of May 26, 1973 and March 20, 
1974 , have been chosen and discuss<'d hero. 

(1) May 26, 1973 

The Doppler record for this event is shown in 
Figure 11 at an operating radio frequency of 4.0125 
MHz. It is evident in this figure that the Doppler 
fluctuations are highly disturbed until 1100 CDT. 

Botii power spectrum density and cross power spectrum 
density analyses indicates that the periodicity of 
the Doppler frequency fluctuations are in the 3 to 5 
minute range. The magnitude of the horizontal phase 
velocities obtained from the cross power spectrum 
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igure 11. CW Doppler record at the operating 
frequency of 4.0125 MHz (event on 
May 26, 1973). 
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density analysis or cross-corrclogram is shown to be 
1000 m/sec propagating: from the south and southwest. 

(2) May 20, 1974 

Figures 12, 13 and 14 are the CW Doppler records 
at operating frequencies 4.0125, 4,759 and 5.734 MHz, 
respectively, under thunderstorm conditions of March 20, 
1974. No tornadoes were sighted in or around North 
Alabama on this day; therefore, it can bo classified 
as a severe thunderstorm active day. The thunderstorm 
associated disturbances are clearly shown oi all the 
records beginning at 1610 CDT and continuing through 
the early morning hours. Wavelike fluctuations with 
long periods of 4 to 5.5 minutes are clearly shown on 
the 4.0125 MHz record beginning at 2200 CDT and lasting 
through 0300 CDT. At higher operation frequencies, 
namely 4.759 and 5.734 MHz, the ionospheric reflected 
signal was lost after 2300 CDT. .The CW Doppler records 
of March 20, 1974, were subjected to cross correlation 
analysis to compute the phase velocities. The horizontal 
trace velocities showed a magnitude of 1000 m/sec 
propagating from a direction southwest of the Doppler 
sounder. lixamination of the Doppler sounder records 
also reveals that the frequency fluctuations have a 
constant amplitude and that they persisted for 4 to 5 
hours . 
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Figure 12. CW Doppler record at the operating frequency of 4.0125 MHz C^vent on March 




c 

o 


c 

(D 

> 


0» 


( 1 

X 

CD 

»o 

r- 

O 

X 

u 

c: 

(U 

3 

cr 

o 

iM 

CA 

fX 


ct5 

<D 

P- 

O 

0^ 

-C 

4^ 


<T$ 

TJ 

• 

O 

U 

o r 

U CD 


5U • 
r-l O 


(X«N| 

a* 

ox 

Q U 

U3?: 


4> 

Im. 

d 

bc 

•r4 

iU 





Figure 14. CW Doppler record at the operating frequency of 
5.743 MHz (event on March 20, 1974). 









B . Tornado Ohservation 

The three tornado events of Novembor 20, 1973, November 
27, 1973 and the extreme tornado outbreak of April 3, 1974 
have been chosen and discussed in this report. 

(1) November 20, 1973 

Figures 15 and 16 show a CW Doppler record for 
November 20, 1973, at the operating frequencies of 
4.0125 MHz, and 4.759 and 5.734 MHz, respectively. 

The wave-like fluctuations are cleat ly seen beginning 
at 1230 CST and continuing tlirough the evening. The 
periodicity of these wave-like disturbances based on 
the analyses from power 'pectrum and cross power 
spectrum is shown to be in the 6 to 8 r.inutc range. 

The phase velocities of the disturbances, derived 
from a cross correlation analysis, arc propagating 
from the south-west with a magnitude of 600-700 m/scc. 
The Doppler record of November 20, 1973, is expanded 
on the time axis to clearly show the time variations 
(Figure 17). It can be seen from Figure 17 that the 
Doppler fluctuations are not quas i - sinuso idal but 
are S-type waves with several Doppler fold backs. Radar 
weather summary charts for Ncvcml'cr 20, 1973 are shown 
in Figure 18. The wcatner map clearly sho^^•s the areas 
of radar echoes, line^ of radar echoes and heights 
of echo tops. The liius of radar echoes arc from the 
southwest to the north of the Doppler system. The 
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Figure 15. CW, Doppler record at the operating frequency of 
4.0125 MHz (event on November 20, 1973). 
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Figure 18. 
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heights of echo tops are of the order of 40,000 to 
50,000 feet (or about 12 to 15 Km). The radar summary 
charts indicate that the severe weather activity is 
located southwest of our Doppler array, which is the 
direction from which the wavelike disturbances seen 
on Doppler records are propagating. 

2. November 27, 1973 

Figure 19 shows the Doppler records at the operating 
frequency of 4.0125 MHz on November 27, 1973 under 
tornado conditions. This figure clearly shows the 
wavclikc fluctuations beginning about 1100 CST and 
continuing through the evening hours. The period of 
the fluctuations obtained from both power spectrum 
and cross power spectrum analysis arc indicated to be 
on the order of 6 to 8 minutes. The horizontal trace 
velocities of the disturbances , derived from a cross 
correlation analysis, are also propagating from the 
southwest with a magnitude of 600-700 m/sec. The 
radar weather summary charts for November 27, 1973, 
are shown in Figure 20. The weather map also clearly 
shows that the lines radar echoes are from the south- 
wes' to the north of the Doppler system. The heights 
of echo tops are of the order of 40,000 to 50,000 feet 
(or about 12 to 15 Km). Thus, the radar summary charts 
indicate that the severe weather activity is located 
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Figure 19. CW Doppler record at the operating frequency of 4.759 
and 4.0125 MHz (event on November 27, 1973). 
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southwest of our Doppler array, which is the direction 
from which the wavelike disturbances seen on Doppler 
records are propagatinp. The oscillations of Doppler 
fluctuations in this case persist more than 6 hours. 

(3) April 3, 1974 

During the extreme tornado outbreak of April 3, 
1974, we observed wave-like fluctuations other than 
the ordinary traveling ionospheric disturbances (TID) . 
Figures 21, 22 and 23 show the CW Doppler records for 
this day at the operating frequencies of 4.0125, 4.759 
and 5.734 MHz, respectively. The observed wave-like 
fluctuations show several S-shaped Doppler traces or 
Doppler frequency "fold backs". The periods of these 
atmospheric waves can he determined by auto-correlating 
the Doppler records. Figures 24, 25, 26 and 27 show 
the auto-spectral density of the v;ave periods observed 
during different time periods. In Fi gure 24 three peaks 
corresponding to waves with periods of 3 min., 6 min., 
and 11 min. arc apparent in the data for the 0700 to 
0800 CDT time period. As shown in Figure 4, at 180 Km 
altitude, the Brunt-Vaisala wave-period is about 10 min. 
Therefore, waves with 3 min. and 6 min. wave periods are 
acoustic waves and waves with 11 min. wave periods are 
gravity waves. The horizontal phase velocity of the 
disturbances, calculated from a cross-correlation 
analysis, shows they are propagating in n direction 






PAGB Q 
POOR QUAIJTjy 


Figure 21. CW Doppler record at the operating frequency 4.0125 
(event on April 3, 1974). 
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Figure 22. CW Doppler record at the operating frequency 4.759 mz (event on April 
1974) . 
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Figure 23. CW Doppler record at the operating frequency 
April 3, 1974) . 
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Figure 24. Auto spectrum density of the ionospheric 

disturbance;-, on April 3, 1974 at 0700-0800 CDT. 
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Figure 25. Auto spectrum density of the ionospheric disturbances 
of April 3, 1974 at 1200*1300 CDT. 
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Figure 26. Auto spc'ctrun density of the ionosphorie disturbance.- 
of A|.ril 3, 1974 at 1300-1400 CDT. 
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103“ from goomaj»nctic north at a speed of 225 m/scc 
which is a typical pliasc velocity for gravity waves at 
F-region ionospheric heights. Gravity waves are 
predominating over infrasonic wave-; during 0700 - 0800 
CDT time period. Possible 'sources of the disturbances 
will be discussed in later sections. 

The horizontal phase velocity of the disturbances 
can be calculated front cross -correlograms . This is 
because the zero frequency crossings are essentially 
independent of the- amplitude of the change of height 
of the reflecting surface and they can be used to 
deternine the average period and average time displace- 
ment, Figure 28 ‘^hows cross-corrcl ograms during the 
time interval 14C0-1.S00 CDT on 3 April 1975. The 
average wave-period shown in this cros‘;-corrclogram 
is 13 min. This is in good agreement with tlie wavc- 
periods dui ing the same time interval obtained from auto- 
spectral analysis, as shown in Figure 27, The time 
displacements of the signals between two stations 
(height peaks from correlograms) arc as follows: 


'''T) F-M 

■= 650 

sec 


(‘T)f.n 

— 330 

sec 

(5-1) 


H 

I ^ 

o 

see 



where the subscript F-M denotes the time delay between 
Ft. McClellan and Musclt- Shoals; subscript F-N, time 
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Delay Time 

Figure 28. Cross correlograms of the rero crossing from N'ickajack 
Dam fN), Muscle Shoals fM) and Ft. McClellan fF) . 


52 



f 


delay between Ft. McClellan and Nickajack Dam; and 
subscript N-M, time delay between Nickajeak Dam and 
Muscle Shoals. Acceptable combinations of time 
displacements must satisfy the relat ion.ship 

~ (5*2) 

In our case, it can bo shown that the highest peaks 
of the corrologram, as shown in equation (S-1), satisfy 
Equation fS-2). The horizontal pliase velocities (both 
magnitude and direction clockwise from geomagnetic 
north) during various time intervals on April 3, 1974 
are shown in Table 3. 

Radar summary charts provided by the National 
Weather Service, Figures 29, 33, 31 and 52 were used 
in comparing results calculated from both auto- 
spectral analysis (Figures 24 to 27) and cross- 
correlograms (Figure 28 and Table 3) show line echoes 
and cloud tops from werither r.adar summaries at di fferent 
time intervals. During the time interval 0700 to 0800 
CDT, Figure 24 and TaHe 3 show that gravity waves 
with wave periods of 11 min. and horizontal phase 
speeds of 225 m/sec propagating from the northwest 
(towards 103*’ from north). Figure 29 suggests that 
the (1 i si url>ances mav have' originaied in the severe 
storm cells in tlu' lines of radar echoes with highest 
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DOPPLfR FLUCTUATIONS OBSERVED DURING SEVERE WEATHER ACTIVITY 
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Figure 29. '^^dar su'”'^ary c^art cf 3. 1^7^ ^T 0735 C^T. 






Figure 30, Radar summary chart of April 3, 197^1 at 1235 CDT. 





Figure 32. Radar summary chart of April 3, 197A at 1535 CDT. 


echo tops from 30,000 to 48,000 ft. (or 10 to 15 Km). 
During the time interval 1200 to 1300 CDT, Figure 25 
and Table 3 indicate that gravity waves with wave 
periods of 12 min. and horizontal phase speeds of 
101 m/sec were prop.agat ing from the northeast (toward 
252® from north). Figure 30 indicates that the source 
of these disturbances may be in the severe storm cells 
which are shown in lines of radar echoes with highest 
echo tops from 47,000 to 50,000 ft. (or 14 to 16 Km). 
During the time interval 1300 to 1400 CDT, Figure 26 
and Table 3 indicate that gravity waves with wave 
periods of 13 min. and horizontal phase speed of 217 
m/sec were propagating from the northeast (toward 
225® from north). Figure 31 shows that the source 
of these disturbances may be in the severe storm cells 
which arc shown in lines of radar echoes with highest 

echo tops from 37,000 to 56,000 ft. (or 11 to 16 Km). 

During the time interval 1400 to 1500 CDT, Figure 27 
and Table 3 show that gravity waves with wave periods 
of 12 min. and horizontal phase speeds of 16C m/sec 
were propagating from the northwest (toward 146® from 
north) . 

Again, Figure 32 indicates tiic source of these 
disturbances may be in the severe storm cells which 
are shown in lines of radar echoes with the highest 

echo tO[)S from 50,000 to 62,000 ft. (or 15 to 19 Km). 
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VI. DISCUSSIONS AND CONCLUSIONS 


Georges (1968; 1973), Baker and Davies (1969), and 
Davies and Jones (1972) reported observation of ionospheric 
disturbances during severe storms and suggested that infrasonic 
waves could be generated by the storm system. During the past 
two years, the experimental observations accomplished by our 
group also indicate that infrasonic waves with wave periods 
of 2 to 8 minutes and with horizontal phase velocity 600- 
1000 m/sec which is in the range of acoustic velocity at 
ionospheric height are generated by severe thunderstorm and 
tornado activity. Furthermore, during the tornado outbreak 
of April 3, 1974, we observed gravity waves with wave periods 
of 11 to 13 min. and horizontal phase speeds of 100 to 225 m/sec 
with several Doppler fold backs in tlie Doppler records. This 
result is slightly different from those of the Boulder group; 
in particular, the gravity waves in this case had wave periods 
of 11 to 13 min. as compared with the infrasonic waves with 
wave periods of 3 to 5 min. in the observations by the Boulder 
Group. However, recent observations made by a Canadian grouo 
at London, Ontario, also contained thunderstorm induced gravity 
waves with wave periods of 12 to 20 min. (Curry and Murry, 
1974). Similarily, French and U.K. groups also reported 
observations of gravity waves with wave periods of the order 
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of 20 min. in the vicinity of a weather disturbance fRcrtin, 
et. al . , 1975) . 

It is possible that waves with periods of 2 to 8 min. 
are transformed to waves with periods of 11 to 13 min. by the 
superimposition of Doppler shifts due to neutral winds by the 
following relation 

w' = w - k • Ij (6-1) 

where w is the wave frequency observed on the ('.round; the 

intrinsic frequency as seen by the observer moving with the 
wind; k, the wave vector of the disturbances; and U, the wind 
velocity. However, waves with periods of 2 to 8 min. at 
ionospheric heights belong to the category of acoustic waves 
while those with periods of 11 to 13 min. are in the category 
of gravity waves. This transformation of acoustic waves to 
gravity waves is possible only when the wave-wave interaction 
dominates the dynamics of up]>er atmosp’ ore fSagdeev and Galeev, 
1969). The super- imposi t ion of a Dopplf-r shift due only to 
neutral winds cannot transform iiifrasonic waves into c.ravity 
waves. This seems to indicate that: (1) severe weather 

activity is responsible for the generation of both gravity 
waves and infrasonic waves; and (2) the dynamics of convective 
motion, wind shear and/or other unknown mechanisms might be 
rcs])onsible for the process of wave-wave inter- act i on required 
to transform acoustic waves which were generated by severe 
weather activity to eravity waves. 
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The waveforms of the Doppler records shown in Figures 15, 
16, 17, 19, 21, 22 and 23 are highly distorted in the S-shaped 
fasion. This S-shaped record has been explained in terms of a 
corrugated reflecting surface that moves horizontally with 
constant speed (Davies and Jones, 1972). We have assumed that 
the transmitter and receiver both move with speed V while the 
ionosphere is stationary. If the amplitude of the corrugated 
surface is A, the mean height of the reflecting surface is H, 
the path, S, from the transmitter to the reflecting point is 
given by (Davies and Jones, 1972) 

1 /2 

S = (H + A cos 9) (1 + k* A^ sin^0) (6-2) 

where 9 = kx is the phase of the reflecting point on the 
corrugated surface, and x is the horizontal spatial axis which 
the transmitter moves along. The relation between phase, 9, 
and time, t, is given by 

^ - kA (H + A cos9) sin 9 ® Vt (6-3) 

Thus, the instantaneous Doppler shift, Af, for a carrier 
frequency, f, is: 


Af _ 
T" “ 


C 


sin 0 


(1 + k* A* sin*9) 




X [1 • k*IIA cos 0 - k*A*cos 20] 


(6-4) 


62 



1 


f 


Davies and .lones (1972) dropped the second and tliird terms 
in the square bracket in liquation (6.4). It will be shown 
later that the second term should have been retained since 
it is on the order of unity which is essentially the same 
order as the first term. 

For small amplitude waves, <'< 1, and Fquation (6-2) 

becomes 


II 2 

cos*0 ♦ ( 7 f) cos*0 - cos 0 + 


1 *a2 

k A 


^ -rS - H) = 0 

k^ll’ 


(6-5) 


which is a third order ilgcbraic equation of cos 0. In Rcneral , 
Hquation (6-5) has three roots which means that three echoes 
at different phases can he observed at the same time. However, 
for a loiiRcr wavelength wave, k is shorter. Thus, 2/k*A* >>>1, 
and liquation (6-5) reduces to a first order algebraic equation 
since |cos e| £ 1. From Kquation (6-5) the condition for the 
existence of multiple-valued echoes is 

k^ AM ^ 1 (6-6) 

When this condition exists, the possibility of observing S-shaped 
Doppler traces is solely dependent on decreasing the wavelength 
and/or increasing the amplitude of the disturbances. 

Hquations (6-2), (6-3) and (6-4) have been solved 
numerically. Figure 33 shows phase path vs. time to reach 
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the receiver . It is clear that the phase path is sinusoidal 
when the ratio of the amplitude of the wavelength is less than 
0.03. When the ratio of the amplitude to the wavelength is 
greater than 0.03, it is possible to observe a three-valued 
phase path. Figure 34 shows the Doppler phase shift, Af, vs. 
time to reach the receiver. When the ratio of the amplitude 
to the wavelength is less than 0.03, the Doppler phase shift 
is practically quasi-sinusoidal . Doppler fold-backs or 
S-shaped Doppler traces arc formed when the ratio of the 
amplitude to the wavelength is greater than 0.03 for this 
case . 

In conclusion, the analysis of the Doppler records 
during the severe weather activity reveals the following 
results: (1) both infrasonic waves with wave-periods of 2 

to 8 min. and with horizontal phase velocity of 600-1000 
m/sec, and gravity-waves with wave-periods of 11 to 13 min. 
and horizonnl phase velocity of lOf' to 225 m/sec arc observed 
during severe thunderstorm activity and tornado activity; (2) 
several Doppler fold-backs are observed during the tornado- 
active day and quasi-sinusoidal wave -prof ilcs arc observed 
during the severc-thunderstorm-active day; f3) the Doppler 
fold-backs phenomenon is due to radio wave reflection from 
disturbances whore the amplitude tn wav(‘lcngth ratio is 
high; (4) the eastward wind hetwetn the source and the observed 
fluctuations might be responsible for the change of the wave- 
period from infrasonic waves to grav i t \ -waves . 
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